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Abstract 

The rates are measured per hadronic decay for gluon splitting to bb quark pairs, g^_^, and 
of events containing two bb quark pairs, (74b, using a sample of four-jet events selected from data 
collected with the OPAL detector. Events with an enhanced signal of gluon splitting to bb quarks 
are selected if two of the jets are close in phase-space and contain detached secondary vertices. 
For the event sample containing two bb quark pairs, three of the four jets are required to have 
a significantly detached secondary vertex. Information from the event topology is combined in a 
likelihood fit to extract the values of and (74b, namely 



g^-^ = (3.07±0.53(stat)±0.97(syst)) X 10"^ 
(0.36 ± 0.17(stat) ± 0.27(syst)) x 10"^. 



34b = 
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1 Introduction 



Bottom quark pairs in TF' decays can be produced either directly via bb or indirectly, when a 

gluon is radiated from a quark and then splits into a bb quark pair. A special case consists of events 
with both direct and indirect b quark production, — > bbg bbbb. The rates and (74b per 
hadronic decay, for the reactions Z° qqg with g — > bb and Z° — > bbbb are sensitive to both the 
b quark mass and the strong coupling constant a^. Hence measurements of these rates are tests of the 
theory of Quantum Chromodynamics (QCD). Events with gluon splitting into bb quark pairs are an 
important background for the measurement of i?b) the fraction of hadronic Z*^ decays into bb quark 
pairs. Therefore, a more precise determination of might lead to a reduction of the systematic 
uncertainty in i?b- 

The rate has been calculated in ||l|, ^ , including the re-summation of leading logarithmic terms. 
Those authors point out that the parton-shower approach as implemented e.g. in Jetset 7.4 [^] is a 
good approximation. Numerical calculations of g^^ are given in ^, Q|, predicting a rate in the range 
^theor _ ^2 g _ 2 X 10~^, depending on the b quark mass and the strong coupling constant. 

An estimate of (74b can be obtained from the rate for direct bb production, i^b, multiplied by the 
rate of indirect bb production, gT^. This simple picture is modified because the phase-space for two bb 
quark pairs is smaller than that for two light and two b quarks. The interference between secondary 
bb production and primary bb production cancels to zero at leading order in cks, except for the bbbb 
final state 0, ^. In |jl|, the contribution of this interference term is shown to be less than 0.2% of g^^. 

Measurements of g^ using four-jet final states have been reported by the DELPHI and ALEPH 
collaborations, with the results = (2.1±l.l(stat)±0.9(syst)) x IQ-^ |] and s-j^^ = (2.77±0.42(stat)± 
0.57(syst)) X 10~^ ||^, respectively, in agreement with theoretical predictions. A measurement of the 
bbbb final state in a three jet topology has recently been presented by DELPHI [P, with the result 
54b = (0.60 ± 0.19(stat) ± 0.14(syst)) x 10~^. This result is translated into a gluon splitting rate 
5bb ~ ^ l.O(stat) lb 0.8(syst)) x 10~^, using a tree-level QCD calculation for the ratio g^^j g^h- 
The calculation was carried out by DELPHI P], using the Wphact 1.3 Monte Carlo generator 

In this analysis, decays of the into four-jet final states are investigated. Jets from b or b 
quarks are identified by reconstructing secondary decay vertices. The invariant mass of bb quark pairs 
originating from gluons tends to peak just above threshold. This leads to a small relative momentum 
of the two b hadrons produced in the fragmentation process. By contrast, directly produced bb 
quark pairs have high invariant masses, since they carry a large part of the Z*^ energy. This and other 
characteristics are used to select event samples enriched in the process g — > bb. An angular correlation 
defined similarly to the Bengtsson-Zerwas angle |10| is used to further differentiate between qqqq and 



qqgg final states. In addition, events with three reconstructed secondary vertices are selected. They 
are used to measure the g'4b rate. 



2 Event selection and reconstruction 



2.1 The OPAL detector 



The OPAL detector is described in detail elsewhere [11 1. Only a brief description of the detector 
elements relevant to this analysis is given here. Charged tracks are reconstructed in the central 
tracking system. It consists of a silicon microvertex detector, a vertex drift chamber equipped with 
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axial and stereo wires, a large jet chamber and z chamber^. A solenoid providing a uniform magnetic 
field of 0.435 T parallel to the z-axis surrounds the central tracking system. The silicon microvertex 
detector [12| has two layers which measure tracks in {r^^\ This detector was upgraded in 1993 to 
provide a precise measurement of the z coordinate [O]. Before this detector was upgraded again in 



1995 for the high energy operation in the LEP 2 programme [14| the inner layer covered the range 
|cos(0)| < 0.83 and the outer layer the range |cos(0)| < 0.77. The vertex chamber extends over 
the range |cos0| < 0.95. The coil is surrounded by scintillators for time-of-flight measurements and 
a barrel lead-glass electromagnetic calorimeter. Including the endcap electromagnetic calorimeter, 
the lead-glass blocks cover the range |cos0[ < 0.98. The magnetic return yoke is instrumented with 
streamer tubes and serves as a hadron calorimeter. The return yoke in turn is surrounded by muon 
chambers. 



2.2 Event selection and reconstruction 

The analysis uses data taken with the OPAL detector in the years 1992-1995 on or near the Z*^ 
resonance. Hadronic decays are selected with an efficiency of 98.4%, as described in |jl5|. Only 
events with the tracking system and the electromagnetic calorimeter fully operational are used in this 
analysis. A total number of 3.35 million hadronic events are selected. In these events well-measured 
charged tracks are used with a momentum pt > 0.15GeV/c in the {x^y) plane, and clusters in the 
electromagnetic calorimeter with energies above 0.1 GeV (0.25 GeV) in the barrel (end-cap) region. 
The energies of clusters pointed to by charged tracks are corrected for double counting by subtracting 
the energy deposition expected from the track momentum |16|. 



2.3 Simulated events 

A total number of 150 000 events generated with Pythia 6.130 |17] and 14 million events generated 



with Jetset 7.4 Q are used to evaluate the efficiencies for signal and background, respectively. 
Within the Jetset 7.4 sample, 4.3 million (1.9 million) events were generated in special runs with 
a primary bb (cc) quark pair to increase the statistical significance of the description of background 
processes with secondary decay vertices. To study the g — > bb signal process 150 000 events with 
gluon splitting to bb in the parton shower were generated with Pythia 6.130. This generator is used 



with a special option for gluon splitting to massive quarks, as explained in Section 4.1. All signal and 



background events were passed through a complete simulation of the OPAL detector |18]. 



The simulated events are weighted to correspond to the measured values of cc and bb production, 
Rc = 0.1724 and = 0.21664, given in [O]. The rate of gluon splitting to cc pairs is set to 
9cc = 3.20 X 10^^ as measured by OPAL pO| . The ratio of the number of events with primary 
produced b quarks and g ^ cc to the total number of events with g — > cc is fixed to the value 
predicted by Jetset 7.4. 

In the following, simulated events with gluon splitting to bb signal are referred to as g ^ bb. 
They are further differentiated into events with a primary charm or light quark, referred to as qqbb 
(q = udsc), and events with two bb quark pairs, referred to as bbbb. 



^The OPAL coordinate system is defined with positive z-axis along the electron beam direction, the x axis pointing 
to the center of the LEP accelerator ring and the y-axis normal to the x-z plane. The polar and azimuthal angles are 
denoted by Q and </!>, respectively. 
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Background to this analysis consists of various sources of four jet events. The dominant four- 
jet process is the production of two gluons in addition to the primary quarks, either from double 
Bremsstrahlung or from the triple gluon vertex. About 7% of the four-jet events are expected to be 
from gluon-splitting to a quark antiquark pair pl|]. 

As this analysis is based on the identification of secondary decay vertices, there are two sources of 
background which have to be studied more carefully: events with gluon splitting to cc accompanied 
by any flavor of primary quark, and four-jet events with a primary bb quark pair but without gluon 
splitting to heavy quarks. They are referred to as bbxx (x = guds). 

For comparisons to the data, the rate of gluon splitting to bb pairs is set to g-^^-^ = 2.5 x 10^^ and 
the rate of bbbb events is set to (74b = 0.4 x 10~^. Note that the rate g^-^ includes bbbb events. The 
fit procedure, described in Section is completely independent of these two numbers. 



2.4 The four-jet selection 



The four- momenta of the selected tracks and clusters are combined to form four jets, using the kj_ 
(Durham) algorithm |2^. The value ycut at which an event makes a transition between a three-jet 
and a four-jet assignment is denoted 2/34. Figure || shows the normalized distribution of the quantity 
2/34 for events that are selected as hadronic decays, thereby comparing the data to the Monte Carlo 
prediction. In addition, the distribution predicted for g ^ bb events is shown, scaled by a factor of 
400. These events populate mainly the region of high 1/34. A cut 1/34 > 0.006 is made to define the 
four-jet sample. This cut rejects nearly 90% of the background events, while retaining about 60% of 
the signal events. The estimated signal fraction assuming a signal rate g^^-^ = 2.5 x 10"'^ is 1.2%. In 
the data 443 334 events are selected, corresponding to 13.0% of all hadronic Z" decays. In the Monte 
Carlo simulation only 11.5% of the events are selected, because the prediction is slightly shifted to 
lower values of 1/34 with respect to the data. This corresponds to a deficit of (11.7 ± 0.2 it 0.8)% 
four-jet events in the simulation compared to the data, where the first uncertainty is due to the finite 
number of events in the data and the simulation, the second is from the uncertainty in the rate of 
gluon splitting to cc, gcc = (3.20 it 0.43) x 10~^ ||2^. This deficit is attributed to missing higher orders 
in the parton shower simulation. To deal with this normalization problem, the number of simulated 
events is normalized to the number of four-jet events throughout this analysis, instead of normalizing 
to the number of hadronic decays. The stability of this analysis with respect to variations of the 
cut in 2/34, affecting this normalization, is discussed in Section 



2.5 Vertex reconstruction 



A primary vertex is reconstructed for each event as described in [^]. Then for each jet a secondary 
vertex is reconstructed, using the full three-dimensional tracking information |24]. All tracks in the jet 



satisfying additional quality cuts on the momentum p > 0.5 GeV/c, the distance to the fitted primary 
vertex in the r/(j) plane do < 0.3 cm and the uncertainty on do of f^do < 0.1 cm are fitted to a common 
vertex. The relatively high momentum cut results in an increased contribution of tracks from the 
decay of b hadrons. Tracks with a large contribution are removed from the vertex and the fit is 
repeated until the contribution of each track is smaller than 4. If three or more tracks remain, the 
secondary vertex is accepted. 

The decay length I and its error a are calculated as the distance of the primary to the secondary 
vertex. The direction is constrained to be parallel to the jet axis. The decay length / is positive if 
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Figure 1: Normalized distributions of the jet resolution parameter 1/34. The data are shown with dots, 
the shape predicted from Monte Carlo simulations with a solid line and the shape predicted for the 
signal events as a dotted line. 

the angle between the jet axis and the vector pointing from the primary to the secondary vertex is 
less than 90 degrees, negative otherwise. Vertices with I > are used to identify b hadrons. In this 
analysis, the following variables are used to identify secondary vertices originating from the decay of 
b hadrons: 

• The decay length significance //o", calculated from the decay length / and the experimental 
uncertainty cr on /. Secondary vertices significantly separated from the primary vertex are 
selected with a cut l/a > 3. 

• For the vertices surviving the l/a > 3 cut the output AW of a neural network is calculated. 
This neural network has been trained to separate vertices originating from b hadron decays from 
those in charm or light quark events. 

The neural network was developed for the OPAL i?b analysis [^]. It has five inputs: the decay length 
significance l/a, the decay length I, the number of tracks in the secondary vertex Ug, the reduced decay 
length Ir/ar, where one well-defined track [^] has been removed from the vertex fit, and a variable 
xd, sensitive to the invariant mass of the tracks in the jet that have a high probability of originating 
from a b hadron decay. The neural network output AW lies between zero and one. Values close to 
one indicate a high probability that the vertex is associated with a b hadron decay. 
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class "2+2" 



class "3+1" 



g — > bb 
candidates 




g ^ bb 
candidates 




4 jets I 3 jets i 2 jets 



2/23 



y34 



2/23 



Figure 2: Illustration of the definition of class "2+2" and class "3+1" and the g ^ bb candidate jet 
selection. 

3 Analysis 



In Section 3.1 the selection of candidate events for gluon splitting is described. By changing the y-cut 
and exploiting the transition from four to three jets, two jets are selected in each event as candidates 
to have originated from gluon splitting. The candidate jets are checked for secondary vertices and 
events with two significant secondary vertices are selected. The event sample is subdivided into two 
distinct classes, depending on the event topology. Optimized cuts on the neural network outputs 



are applied for each class, to define the candidate events. In Section 3^ a dedicated selection of 
candidates for the process —>■ bbbb is discussed, where all four jets are checked for secondary 



vertices. Finally, in Section the rate of gluon splitting to bb is calculated. For each of the selected 
event samples, angular distributions sensitive to four-quark final states are studied. The rates g^^-^ 
and 5'4b are calculated using a binned maximum likelihood fit. The signal and background selection 
efficiencies and these angular distributions are used as input to the fit. 



3.1 The qqbb event selection 

In each four-jet event the y-cut is increased, until the event changes to a three-jet event (ycut > 2/34)- 
The two jets that are combined in this step are considered as candidates for g bb. The y-cut is 
then increased further, until the event changes to a two-jet event (i/cut > 2/23)- There are two distinct 
possibilities for this, as shown in Figure 

Class "2-|-2": Events belong to this class if none of the original four jets is identical to any of the 
two jets obtained after increasing the y-cut to force the event into a two-jet event. 

Class "3-|-l": Events belong to this class if one of the original four jets is identical to one of the two 
jets obtained after increasing the y-cut to force the event into a two-jet event. 
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Number of events 


Selection 








sig 


nal 


efficiencies 




data 


bgnd 


qqbb 


bbbb 


^qqbb 


''bbbb 


The four-jet selection 


Total four-jet 


443334 


438191 


4263 


880 


52.72% 


56.95% 


"2-h2" 


235201 


238500 


1901 


373 


23.51% 


24.15% 


"3-hl" 


208133 


199691 


2362 


507 


29.21% 


32.80% 


Event selection in Class "2-1-2" 


(V^)2 > 3 


613 


596.2 


51.2 


17.8 


0.63% 


1.15% 


NNi + NN2 > 1.7 


39 


25.8 


8.6 


3.0 


0.11% 


0.20% 


Sample A 


14 


6.3 


1.2 


2.2 


0.02% 


0.15% 


Sample B 


25 


19.5 


7.4 


0.8 


0.09% 


0.05% 


Event selection in Class "3-1-1" 


{1/(7)2 > 3 


359 


310.6 


71.4 


22.3 


0.88% 


1.31% 


NN1 + NN2 > 1.1 


153 


92.8 


42.9 


13.2 


0.53% 


0.85% 


Sample C 


44 


40.1 


6.7 


9.6 


0.08% 


0.62% 


Sample D 


109 


52.7 


36.2 


3.6 


0.45% 


0.23% 


The dedicated bbbb selection 


{l/(Th > 3 


628 


642.8 


18.0 


55.8 


0.22% 


3.61% 


remove overlap with A-D 


589 


604.6 


14.0 


46.0 


0.17% 


2.97% 


Sample E 


29 


21.1 


0.3 


9.1 


0.004% 


0.59% 


Selected (A E) 


221 


139.7 


51.8 


25.3 


0.64% 


1.64% 



Table 1: Events selected at the different steps of the analysis (data), and the number of background 
(bgnd) and signal events expected from the simulation with g^-^ = 2.5 x 10^'^ and (74b = 0.4 x 10^'^. 
Also shown are the efficiencies to select qqbb (q = udsc) and bbbb signal events, predicted from the 
Monte Carlo simulation. 

Monte Carlo studies show that events from the process g ^ bb are preferentially selected in class 
"3-|-l" . The corresponding selection efficiencies are given in Table For a signal rate ^^j^ = 2.5 x 10^^, 
in class "3-1-1" 1.47% of the events are from g — > bb. In class "2-1-2" only 0.99% of the events are 
from g —5- bb. The event classification procedure selects 235 201 data events in class "2-1-2" and 
208133 in class "3-1-1". The Monte Carlo simulation predicts 240 774 and 202 560 events, respectively. 
Implications on the final result for g^^ and (74b from this statistically significant difference in the 
number of observed to expected events in those two event classes are considered in the discussion of 
the experimental uncertainties. 

The two g — 5- bb candidate jets selected in the first step are checked for secondary vertices. If both 
of these jets have a reconstructed secondary vertex with l/a > 3, the event is selected. This cut is 
studied in Fig. |3|a, where the decay length significance {l/cr)2 of the second g — > bb candidate jet is 
shown. The two jets are ordered such that {l/cr)i > (Z/cr)2. If the variable {1/ct)2 is larger than 3, the 
event is selected, otherwise it is rejected. Note that {l/o')2 > 3 implies {l/o')i > 3. Because only the 
smaller decay length significance of the two vertices enters this distribution, there are many entries at 
negative {l/cr)2 values. The cut {l/a)2 > 3 reduces the background from events with light flavors and 
it also reduces the fraction of events from the g — > cc process. 

In the next step, the output of a neural network trained to recognize vertices from b hadron 
decays is calculated for the reconstructed vertices in the selected events. The two g ^ bb candidates 
are ordered by their corresponding neural network output AW such that NNi > NN2- The sum 
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Figure 3: Variables used to select g — > bb candidate events. Plot (a) shows the decay length signif- 
icance {l/cr)2 of the g — > bb candidate jet with the smallest decay length significance. Plot (b) and 
(c) shows the sum of the neural network outputs NNi and NN2 for the events surviving the l/a cut 
in class "2+2" and "3+1", respectively. The cut values are indicated by arrows. The points show 
data and the histogram shows the Monte Carlo simulation. The contributions from various four-jet 
processes are indicated. 



NNi + A'7V2 shows good sensitivity to the signal process g — > bb for both event classes "2+2" and 
"3+1". It is depicted in figures |b and ^. To enrich the signal, cuts NNi + NN2 > 1.7 in class 
"2+2" and NNi + NN2 > 1-1 in class "3+1" are made. The cut values are chosen such as to obtain 
purities of about 40% for each class of events. Some discrepancies between data and the prediction 
are observed in Fig. ^ and in table |l|, mainly at low values of l/a and low values of NNi + NN2. These 



can be explained by uncertainties in the knowledge of the detector resolution, studied in Section 4.3. 



Finally, the other two jets that were not considered to be gluon splitting candidates are checked for 
secondary vertices to gain some separation power for the bbbb signal events. These jets are referred 
to as "primary quark jet candidates". If both of these jets have a vertex, only the jet with the larger 
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Figure 4: Largest decay length significance {I/ct)^ of the primary quark jet candidates, as defined in 
the text, after applying all selection cuts for the qqbb selection. Events from class "2+2" appear in 
Fig. 0a, events from class "3+1" in Fig. ^b. The points show data and the histogram shows the Monte 
Carlo simulation. The contributions of various four- jet processes are indicated. The cut at {l/o')p = 2 
separating the event samples A from B and C from D is indicated. 



decay length significance is considered. Fig. ^ shows the distribution of the decay length significance 
of this reconstructed vertex, {l/a)p, for the events selected from class "2+2" and "3+1". The variable 
(Z/(t)p shows some separation between the bbbb events and the qqbb (q = udsc) events. A cut at 
(//(t)p = 2 is made to select event samples enriched or depleted with bbbb events. The following four 
event samples are defined 



Event sample A: Events from class "2+2" with NNi + NN2 > 1.7 and {1/(t)p > 2, enriched in 
bbbb. 

Event sample B: Events from class "2+2" with NN1+NN2 > 1.7 not selected in sample A, depleted 
in bbbb. 

Event sample C: Events from class "3+1" with NNi + NN2 > 0.7 and {l/cr)p > 2, enriched in 
bbbb. 

Event sample D: Events from class "3+1" with NN1+NN2 > 0.7 not selected in sample C, depleted 
in bbbb. 
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□ qqxx (q=udsc,x=guds) 



-f- DATA ^ bbxx (x=guds) 

— MC g-^ cc (primary udscb quark) 

^ bb (primary udsc quark) 
m g^ bb (primary b quark) 

Figure 5: Variables used to select bbbb events. Plot (a) shows the third largest decay length sig- 
nificance {l/(T)-i of the four jets, plot (h) the third largest neural network output NN^ with the cut 
(V^)3 > 3 applied and the overlap with the event samples A-D removed. The points show data and 
the histogram shows the Monte Carlo simulation. The contributions from various four-jet processes 
are indicated. Arrows indicate the position of the cuts. 

3.2 The dedicated bbbb event selection 

Starting again with the entire sample of four-jet events, selected with y^^ > 0.006, a dedicated 
selection of bbbb events is set up, independent of the cuts presented in the previous section. To find 
events with four b hadrons, events are selected where at least three decay vertices with a decay length 
significance l/a > 3 are found. Fig. |5|a illustrates this selection of three significant vertices, showing 
the third largest decay length significance, denoted {l/a)^. The bbbb signal is enhanced for large 
values of this variable. The cut {l/cr)^ > 3 suppresses light flavors, most of the g ^ cc events, and also 
the qqbb (q = udsc) events. Fig. ^b shows the third largest neural network output AW3 for all selected 
vertices, where the events already selected in sample A-D are excluded to avoid double counting. The 
background dominates the region of low while the bbbb signal extends to high values of NN3. 

A cut NN3 > 0.7 is chosen to select the final bbbb candidates from distribution, denoted as sample 
E. 

3.3 Calculation of g^-^ and gr4t, 

After applying all cuts, 221 events remain in the event samples A-E, where the simulation with 
g^^ = 2.5 X 10~^ and g^y, = 0.4 x 10~^ predicts 217 it 7 events. The efficiencies for selecting signal 
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events, qqbb (q = udsc) or bbbb events, are denoted and e^bbb' '^^^ total efficiencies as obtained 
from the simulation are Cqq^E ~ (0-64 ± 0.02)% and e^bbb ~ (1-64 ± 0.07)%. The uncertainties quoted 
on the number of expected events and the selection efficiencies are due to the limited number of Monte 
Carlo events. Table |l] summarizes the number of selected events for the cuts applied, the number of 
signal and background events expected from the simulation, and the efficiencies to select signal or 
background reactions. 

For the 221 signal events the angle 012-34 between two jet-jet planes is studied to distinguish signal 
events from the background, dominated by events with two quarks and two gluons. The first plane is 
spanned by the two jets that are joined into one jet by the jet-algorithm at the transition from four 
to three jets. The other plane is formed by the other two jets. The definition of this angle 012-34 
is similar to the angular correlation proposed in |l^ to measure the QCD color factors. In Fig. ^, 
the ai2-34 distribution is shown to be consistent with the theoretical prediction. The signal appears 
preferentially at high values of 012-34 while the background has a flatter distribution. 

A maximum likelihood fit of the 26 bins in 012-34 shown in Fig. ^ is performed to extract g^^^ and 
54b, assuming Poisson distributions calculated from the signal and background efficiencies for each 
bin. The likelihood function is given by 

26 

— lnC = ^(/ij — di In/Uj) + const, 

i=l 

where di is the number of observed events in bin i and fii is the corresponding number of predicted 
events (assumed to be the mean of the Poisson distributions). The latter is given by 

(1 - 9cc - gbb)4q + gcc6^qcc + (ffbb " 9ih)ei-^-^; + g4b6tbbb 
fJ-i ^^4-jet 4_jet 4-iet , / n 4-jet , 4-iet ' 

(1 - 5cc - 5bb)eqq + 5cceqq4 + (^bb " 54b)eq-^^ + 94he^-^^-^ 

where A^4-jet is the total number of selected four-jet events used for this analysis. The rate of gluon 
splitting to cc per hadronic decay is denoted ^cc- The two rates g^^-^ and 5f4b are taken as the fit 
parameters. The probabilities to select the signal process for primary light quarks and the bbbb events 
in bin i are denoted e*_, t- and e*T-, t-, while el,-:^-, and e*-; are the probabilities to select the g ^ cc 

qqbb bbbb qqcc qq -t o 

events or other background events in bin i. Finally e'^"'^^*, ef^,'^?-*, e^^j^*, et^^^^ are the efficiencies for 

qqbb bbbb 44'-'' 44 

events in the signal and the background channels to pass the four-jet selection. The likelihood analysis 
leads to the result 

g^-^ = (3.07 ± 0.53) X 10-3, 
54b = (0.36 ± 0.17) X 10-3, 

with a correlation coefficient between g^-^ and (74b of +0.007. The errors are statistical only. 



4 Evaluation of systematic uncertainties 

Various sources of systematic uncertainties are considered, as summarized in Table ^ and discussed in 
the following. Additional cross-checks are presented in Section ||. 



4.1 Model dependence 

Though the parton shower approach as implemented in Jetset 7.4 is expected to describe the gluon 
splitting process quite well |jl], it is desirable to look at various alternative models and at exact 
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Figure 6: The angle 0:12-34 between the plane formed by the two g — > bb candidate jets and the plane 
formed by the other two jets for the g — > bb candidate events after applying all cuts. Figures a-e 
correspond to the event samples A~E. The data are shown as full points with error bars, the Monte 
Carlo simulation as a solid line. The contributions from various four-jet processes are indicated. 



calculations. In Pythia 6.130 [17], using the option MSTJ(42)=3, the calculation of the opening 
angle in the gluon splitting process has been modified to take mass effects into account, as compared 
to Jetset 7.4. The Pythia 6.130 prediction for the g ^ bb process is used to evaluate the main 
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Source of systematic error 




A54b X 10^ 


Correlation 


Model dependence 


0.72 


0.17 


+1 


b quark mass 


0.23 


0.03 


+1 


Monte Carlo statistics 


0.28 


0.09 


-0.07 


Detector simulation 


0.23 


0.15 


+1 


Event classification 


0.18 


0.03 


+1 


Flavor composition Rt~^'^^ 


0.07 


0.02 


+1 


Gluon splittina; to charm Qcc 


0.05 


0.01 


+1 


Bottom fragmentation 


0.39 


0.10 


+0.86 


Bottom decay multiplicities 


0.11 


0.03 


+1 


Bottom production rates 


0.11 


0.01 


+0.70 


Bottom hadron lifetimes 


0.11 


0.05 


+0.09 


Charm fragmentation 


0.07 


0.03 


-0.52 


Charm decay multiplicities 


0.03 


0.01 


+0.96 


Charm production rates 


0.02 


0.01 


-0.97 


Charm hadron lifetimes 


0.05 


0.01 


+0.80 


Total systematic error 


0.97 


0.27 


+0.78 



Table 2: Summary of the systematic uncertainties on the rates g-^-^ and g/i\, and the corresponding 
correlations 



results of this analysis. The Herwig 5.91 [25| Monte Carlo generatoijj provides an alternative model 
for the parton shower. Another approach is given by the Color Dipole Model, using the event generator 
Ariadne 4.08 p|. 



The program Wphact implements e"'"e annihilation to four fermions, where the masses of b 
quarks are taken into account for the calculation of the matrix elements. It can be used to study bbbb 



and bbqq final states. Version 1.3 of Wphact^ |27] allows for the possibility to switch off Feynman 
graphs where the b quarks couple to the or 7, enabling studies of the gluon splitting process at 
tree level. For the case of primary light quarks, calculations are available including next-to-leading 
order logarithmic terms [||, ^. The shapes of the Wphact 1.3 predictions are in agreement with 
these calculations in the regions where such a comparison is possible. 

Fig. |7| shows differential distributions of kinematic variables of the gluon that splits to bb, calcu- 
lated with Pythia 6.130,Jetset 7.4, Herwig 5.91, Ariadne 4.08 and Wphact 1.3. The "gluon" 
variables have been calculated from the bb quark pair at the end of the shower. For the case of bbbb 
events the bb quark pair with the lowest invariant mass is chosen. The hardest energy spectrum is 
predicted by Jetset 7.4, while Ariadne 4.08 leads to the softest energy spectrum, when comparing 
the five models. For the gluon virtuality the Ariadne 4.08 prediction leads to the hardest distribu- 
tion, while Herwig 5.91 predicts the softest spectrum. For the decay angle in the gluon rest frame 
the models differ most significantly at high |cos0*|. The extreme cases are covered by Jetset 7.4 
with a low number of events in this region and by Ariadne 4.08 showing an enhancement of events 
with increasing | cos9*\. Note that the Pythia 6.130 prediction is always well between those extreme 
cases, this is why we decided to use it for the main analysis results. 

To study hadronisation and detector effects, we use two sets of events generated either with 
Jetset 7.4 or with Pythia 6.130, and having a detailed hadronisation and detector simulation. 
The Jetset 7.4 sample is only used for systematic checks. 

^This version was already used in Q. 
^Provided to us by courtesy of the author. 
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Figure 7: Normalized distributions of three kinematic variables of the gluon that splits to bb for 
different models and calculations, as described in the text: (a) the energy E , (b) the virtuality m, (c) 
the decay angle \ cos 9*\ of the bb pair in the gluon rest frame. 



The differential efficiencies to select the signal events, qqbb (q = udsc) and bbbb in this analysis, 
are shown in figure ^. They are evaluated with the Pythia 6.130 event sample. The efficiencies are 
low for small gluon energies E and gluon virtualities m, as it is difficult to detect the two b hadrons 
in two separate jets in this region. The same is true for decay angles | cos9*\ close to one, because in 
this case the two b hadrons have a small transverse momentum relative to the gluon flight direction. 

To evaluate model-dependent effects, the complete analysis chain is performed using either the 
Pythia 6.130 or the Jetset 7.4 event sets. In addition, the Pythia 6.130 events are weighted such 
as to reproduce the gluon distributions from Fig. |^ predicted by Herwig 5.91, Ariadne 4.08 and 
Wphact 1.3, and the analysis is repeated. Table ^ summarizes the individual results of these studies. 
To study possible deficiencies in the reweighting procedure, it is repeated based on the Jetset 7.4 
event sample rather than on the Pythia 6.130 event sample. This results in slightly lower values 
of g^^ and larger values of 5(4b, as compared to the Herwig 5.91, Ariadne 4.08 and Wphagt 1.3 
values shown in table |3[ The most significant differences to the standard reweighting procedure are 
observed for Herwig 5.91, namely g^^^ = 2.88 x 10~^ and g^^ = 0.53 x 10~^. 
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Figure 8: Selection efficiencies as a function of kinematic variables of the gluon that splits to bb; (a) 
the energy E, (h) the virtuality m, (c) the decay angle \ cos0*| of the bb pair in the gluon rest- frame. 
The efficiencies to select qqbb (q = udsc) and bbbb events are indicated. 



Model 


Implementation 


5hh X 


54b X 10^ 


Jetset 7.4 

Pythia 6.130 (default) 
Wphact 1.3 
Herwig 5.91 

Ariad.xe i.08 


Full detector simulation 
Full detector simulation 
Event reweighting 
Event reweighting 

Eveul rowcugliliug 


2.35 ± 0.41 

3.07 ±0.53 
3.22 ±0.58 
3.38 ±0.59 
3. (32 ± 0.6 i 


0.34 ±0.17 

0.36 ±0.17 
0.31 ±0.15 
0.40 ±0.19 
0.36 ±0.19 



Table 3: Results obtained for various Monte Carlo models. 



Finally, the largest differences on g^^ and g'4b to the central value are taken as systematic uncer- 
tainty due to limited knowledge of the gluon splitting mechanism. For g-^^ this is the Pythia 6.130- 
Jetset 7.4 difference, while for g^)^ this is the Pythia 6.130-Herwig 5.91 difference, where the 
Herwig 5.91 result is calculated with the help of reweighted Jetset 7.4 events. 



16 



4.2 Dependence on the b quark mass 



The influence of the b quark mass assumed for the Monte Carlo simulation has been studied 
by changing the b quark mass in the Pythia 6.130 program^ and investigating the properties of 
the b hadrons produced in the fragmentation process. The central value for the b quark mass in 
Pythia 6.130 used throughout this analysis is 5GeV/c^. The shapes of the distributions of various 
kinematic variables are used to reweight the g — > bb Monte Carlo events corresponding to a b quark 
mass of 4.5GeV/c^ and 5.25GeV/c^. This variation covers the uncertainty of the b quark pole mass 
] up to the B meson mass. Note that the way the quark mass parameter is used in Pythia 6.130 
corresponds rather to a constituent quark mass than to a pole mass definition of the quark mass [p9| ], 
which justifies the choice of 5GeV/c^ for the central value. 

The quantities used for the reweighting process are the momenta of the two b hadrons produced in 
the fragmentation process and their invariant mass. These variables are chosen because the efficiency 
to identify b hadrons strongly depends on their momentum. The efficiency to resolve the two b hadrons 
in different jets depends on their invariant mass. The larger deviations from the standard results are 
found for a b quark mass of 4.5 GeV/c^. They are taken into account as systematic uncertainties. 



4.3 Experimental sources of systematic uncertainties 

The modeling of the OPAL detector is important because the analysis depends on an accurate un- 
derstanding of the decay vertex reconstruction. The b tagging efficiencies are mainly sensitive to the 



parameters modeling the production and decay of b hadrons. This will be discussed in Section 4.4. 
The light quark tagging efficiencies in the Monte Carlo simulation are mainly sensitive to details in 
the modeling of the tracking system. Studies are done by increasing the difference of the reconstructed 
track parameters with respect to the true track parameters in the Monte Carlo simulation by 10% 
to cover uncertainties in the knowledge of the flavor tagging variables [^] . This smearing is applied 
separately for parameters defined in the (r, (p) and the (r, z) plane. In (r, (p) a simultaneous smearing 
of the distance of closest approach, do, and the azimuthal angle (po is performed. In (r, z) smearing 
is done simultaneously for the z coordinate of the point of closest approach in (r, (p) and the polar 
angle 9. Finally the gluon splitting analysis is repeated using the modified Monte Carlo sets. The 
uncertainties from the smearing in (r, (p) and (r, z) are added quadratically. 

The normalization to the number of four-jet events revealed differences in the population of the 
event classes "2+2" and "3+1" when comparing the data with the Monte Carlo prediction, as discussed 



in Section |3.1| . This is addressed by repeating the likelihood fit, using the numbers of events and 
selection efficiencies of the "2+2" and "3+1" classes, rather than the number of four-jet events with 
their corresponding selection efficiencies. The variation found is assigned as systematic uncertainty 
due to the event classification. 



4.4 Uncertainties from heavy flavor physics 

The rate of primary bb production in the four-jet sample R^~^^^ can be measured using a double-tag 
technique. Agreement with the Monte Carlo simulation within 2.4% is found, where the statistical 
uncertainties are in the order of 1.5%. The rate R^~^'^^ in the Monte Carlo sample has been changed by 
2.4% and the differences in g^-^ and g^^ to our standard results are taken as systematic uncertainties. 

*The relevant parameter in the Pythia 6.130 program is PMAS(5,1). 
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Varying the rate of primary cc production Rc within its uncertainties only has a negUgible effect on 
the results. 



The rate of secondary cc quark pairs from gluons, Qcc, has been measured by OPAL to be gcc = 
(3.20±0.43) X 10^^ [20|. This measurement is based on lepton identification and the reconstruction of 



D* mesons in three-jet events. It can be considered as statistically independent of the g^^^ measurement 
presented here, which is based on lifetime-tags in four-jet events. The rate gcc thus is varied within its 
uncertainty to obtain systematic uncertainties on g^^ and g^b. The explicit dependence of the results 
of this analysis on ^^cc is given by 

^ = -0.125 x^, ^ = -0.259 x^. 

^bb 9cc 54b 9cc 



The fragmentation functions for charm and bottom quarks are varied to reflect the uncertainties 
in the knowledge of the average scaled energies for D and B mesons, {xe)c = 0.484 it 0.008 and 
{xe)y, — 0.702 lb 0.008 pl|. This is done by varying the parameter e in the parameterization of the 



fragmentation function suggested by Peterson et. al. |32]. In addition the sensitivity to the shape 
of the fragmentation function is checked, using the models suggested by Collins and Spiller |33| and 
Kartvelishvili |3^ . The parameters of these models are chosen to reproduce the measured values of the 
mean scaled energies of the D and B mesons in the Monte Carlo simulation, and the larger deviations 
in g^ and are used as systematic uncertainties from this source. The uncertainties from the 
knowledge of the mean scaled energies and the shape of the distributions are added in quadrature. 



The charged decay multiplicities of the D mesons are varied within the errors given in [p5| around 
the Jetset 7.4 prediction. Particles from decay-chains of b hadrons are excluded from this variation. 
The mean charged decay multiplicity of weakly decaying b hadrons is varied within ub = 4.995 ±0.062 
[^]. This multiplicity includes secondary decays of charm hadrons produced in the decay chains. 
Variations of the neutral decay multiplicities have not been studied, but are expected to be small. 

The production rates of charm and bottom hadrons in the Monte Carlo simulation are varied 



within the ranges given by the LEP Electroweak Working group [31| and the Particle Data Group 
[ p^ ] , respectively. This variation includes particles from primary bb and cc production as well as those 
from gluon splitting to heavy quarks. 

The lifetimes of charm and bottom hadrons are varied around their central value according to the 



numbers given by the Particle Data Group [ 28 1 . 



4.5 Results 



All contributions to the systematic uncertainties are added in quadrature, leading to the results 

g^- = (3.07 ± 0.53(stat.) ± 0.97(syst)) x 10"^ 
54b = (0.36 ± 0.17(stat.) ± 0.27(syst)) x 10~^ 

The systematic errors on g^^ and g(4b have a correlation +0.78. This large positive correlation can be 
understood from the fact that most of the systematic variations influence the selection efficiencies for 
both the g ^ bb and the bbbb signal in the same direction. 
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Figure 9: Dependence of the four-jet rate R^-^et and the fit results on 2/34™, the lower cut in the jet 
resolution parameter 2/34. Figure (a) shows the ratio of the number of events in the data by the number 
of events predicted from the Monte Carlo simulation, figures (h ) and (c) show the fit result for g^^ and 
(74b. The error bars in (b) and (c) correspond to the statistical errors independent from the central 
value yfi'' = 0.006. 

5 Additional cross-checks 



The four-jet selection is defined using a cut 1/34 > 2/34™. Fig. |9|a shows the number of events in the 
data divided by the predicted number of events as a function of this cut. The disagreement in the 



rate of four-jet events, as discussed in Section is clearly visible. The fit results for g^^ and (74b as a 
function of this cut are studied in Fig. ^3 and ^ The fit results are stable under the variation of 2/34™ 
within the independent statistical error. Note that the variation of 2/34™ from to 0.014 corresponds 
to a variation in the number of candidate events from 1003 to 93 and a variation of the estimated 
g — 5- bb signal purity from 10% to 50%. These checks show that the treatment of the normalization 
to the four-jet rate in the likelihood fit is correct. 



To check the fitting procedure a likelihood fit is applied to the normalized signal and background 
shapes in 012-34, calculating directly the fractions of qqbb (q = udsc) and bbbb events in the final 
event selection from these shapes. Using the number of selected events, the total number of hadronic 
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TP decays and the signal efficiencies, and g^h can be calculated from these event fractions. This 
procedure avoids the need to know the absolute selection efficiency for the background. Alternatively, 
the likelihood fit is repeated using only the individual signal and background selection efficiencies of 
the event samples A-E. This analysis is thus independent of the variable q;i2-34- Both results are 
compatible with the default method. 

The fc_L algorithm is used in this analysis not only to define the jets, but also to define the event 
topology. Therefore it is important to test the analysis, using a different algorithm. The JADE EO 
pO| ] jet clustering scheme is used as an alternative method. Four-jet events are selected with a cut 
y^^ADE ^ Q Q]^5^ resulting in a similar number of four-jet events as compared to the default analysis 
which is based on the k± algorithm. To classify the events and define the gluon splitting candidates. 



as described in section |3.l| , and to define the angle ai2-34, as described in Section p.3| , the new 
definition of the jet resolution parameter in the JADE EO scheme is used. All other analysis cuts 
are unchanged with respect to the standard analysis. This results in 207 selected events. The result 
of the likehhood fit using the JADE EO algorithm is g^^^^^ = (3.11 ± 0.59) x lO'^ and gl^^^ = 
(0.45 ± 0.17) X 10^^ (statistical uncertainties only). These results are compatible with the results 
obtained with the k±_ algorithm within the statistical errors, considering the small fraction of events 
common to both selections. 



6 Summary 

A measurement of the inclusive rate of gluon splitting to bb per hadronic TP decay has been performed 
using data taken by OPAL. The result is 

g^- = (3.07 ± 0.53 ± 0.97) x 10"^. 

The rate of events with four bottom quarks is measured simultaneously, with the result 

54b = (0.36 ± 0.17 ± 0.27) X 10"^ 

The correlation of g^-^ and (74b is -1-0.007 from statistical uncertainties and -1-0.78 from systematic 
uncertainties. The ratio 5'4b/9bb ^^'^^ measured to be 

— = 0.116 ± 0.060 ± 0.065. 

fbb 

This is in agreement with the simplified expectation S'4b/5'bb ~ and with a more detailed calculation 
performed by DELPHI |^], predicting g^T^/ g^h = 0.1833 ± 0.0003. The results for g,^ and 54b are 
compatible with previous measurements [y, 0, ^ as well as theoretical predictions [Q, ^ . 
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